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The phase diagram of three-flavor quark matter under compact star constraints is investigated 
within a Nambu-Jona-Lasinio model. Local color and electric charge neutrality is imposed for 
/3-equilibrated superconducting quark matter. The constituent quark masses and the diquark con¬ 
densates are determined selfconsistently in the plane of temperature and quark chemical potential. 

Both strong and intermediate diquark coupling strengths are considered. We show that in both 
cases, gapless superconducting phases do not occur at temperatures relevant for compact star evolu¬ 
tion, i.e., below T ^ 50 MeV. The stability and structure of isothermal quark star configurations are 
evaluated. For intermediate coupling, quark stars are composed of a mixed phase of normal (NQ) 
and two-flavor superconducting (2SC) quark matter up to a maximum mass of 1.21 M©. At higher 
central densities, a phase transition to the three-flavor color flavor locked (CFL) phase occurs and 
the configurations become unstable. For the strong diquark coupling we And stable stars in the 2SC 
phase, with masses up to 1.326 M©. A second family of more compact configurations (twins) with 
a CFL quark matter core and a 2SC shell is also found to be stable. The twins have masses in the 
range 1.301...1.326 M©. We consider also hot isothermal configurations at temperature T = 40 MeV. 

When the hot maximum mass configuration cools down, due to emission of photons and neutrinos, 
a mass defect of 0.1 M© occurs and two final state configurations are possible. 

PACS numbers: 12.38.Mh, 24.85.+p, 26.60.+C, 97.60.-s 


I. INTRODUCTION 

Theoretical investigations of the QCD phase diagram 
at high densities have recently gained momentum due to 
results of non-perturbative low-energy QCD models Ju 
13, 01 of color superconductivity in quark matter 0, 0|. 
These models predict that the diquark pairing conden¬ 
sates are of the order of 100 MeV and a remarkably rich 
phase structure has been identified 0,0 0 0| • The main 
motivation for studying the low-temperature domain of 
the QCD phase diagram is its possible relevance for the 
physics of compact stars [ifl, llJ, [i3| ■ Observable effects 
of color superconducting phases in c omp act stars is ex¬ 
pected, e.g., in the cooling behaviour |l.4ri4ri^ll(Tirr^ . 
magnetic field evolution [la 113 , 12 ( 11 ^, and in burst-type 

phenomena [HEHillir 

The most prominent color superconducting phases 
with large diquark pairing gaps are the two-flavor scalar 
diquark condensate (2SC) and the color-flavor locking 
(CFL) condensate. The latter requires approximate 
SU(3) flavor symmetry and occurs therefore only at 
rather large quark chemical potentials, p-g > 430 — 500 
MeV, of the order of the dynamically generated strange 
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quark mass Ms, whereas the 2SC phase can appear al¬ 
ready at the chiral restoration transition for p-g > 330 — 
350 MeV (2^,113,Ell- Note that the quark chemical po¬ 
tential in the center of a typical compact star is expected 
to not exceed a value of ^ 500 MeV so the volume frac¬ 
tion of a strange quark matter phase will be insufficient 
to entail observable consequences. However, when the 
strange quark mass is considered not dynamically, but 
as a free parameter independent of the thermodynamical 
conditions, it has been shown that for not too large Ms 

the CFL phase dominates over the 2SC phase 113- 
Studies of the QCD phase diagram with fixed strange 
quark mass have recently been extended to the discussion 
of gapless CFL (gCFL) phases mmiii The gapless 
phases occur when the asymmetry between Fermi levels 
of different flavors is large enough to allow for zero en¬ 
ergy excitations while a finite pairing gap exists. They 
have been found first for the 2SC pha se (g2SC) within a 
dynamical chiral quark model jH, ■ 

Any scenario for compact star evolution that is based 
on the occurence of quark matter relies on the assump¬ 
tions about the properties of this phase. It is therefore 
of prior importance to obtain a phase diagram of three- 
flavor quark matter under compact star constraints with 
selfconsistently determined dynamical quark masses. In 
the present paper we will employ the Nambu-Jona- 
Lasinio (NJL) model to delineate the different quark mat¬ 
ter phases in the plane of temperature and chemical po¬ 
tential. We also address the question whether CFL quark 
matter and gapless phases are likely to play a role in com¬ 
pact star interiors. 
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II. MODEL 

In this paper, we consider an NJL model with 
quark-antiquark interactions in the color singlet 
scalar/pseudoscalar channel, and quark-quark inter¬ 
actions in the scalar color antitriplet channel. We 
neglect the less attractive interaction channels, e.g., 
the isospin-singlet channel, which could allow for weak 
spin-1 condensates. Such condensates allow for gapless 
excitations at low temperatures and could be important 
for the cooling behaviour of compact stars. However, the 
coupling strengths in these channels are poorly known 
and we therefore neglect them here. The Lagrangian 
density is is given by 

8 

+ [{qr^qf + {qi^BTfqf] 

a—0 

+ Gd E \{qioif-ijkf-al3^qjfj){qi'a'^i' P'^qj' P') 

fc,7 

j'13')\ ; (l) 

where M°- = diag(TO°, to°, m°) is the current quark mass 
matrix in flavor space and Pij,ap is the chemical potential 
matrix in color and flavor space. Due to strong and weak 
interactions, the various chemical potentials are not in¬ 
dependent. In the superconduc ting phases a [/(I) gauge 
symmetry remains unbroken [ 2 ^, and the associated 
charge is a linear combination of the electric charge, Q, 
and two orthogonal generators of the unbroken SU{2)c 
symmetry. Hence, there are in total four independent 
chemical potentials 

Pij,aP — QpQ)^aP “t” (Ts/is -j- T^p^')6ij ^ ( 2 ) 

where Q = diag(2/3, —1/3, —1/3) is the electric charge 
in flavor space, and T 3 = diag(l,—1,0) and Tg = 
diag(l/v^, 1/V3,- 2 /v^) are the generators in color 
space. The quark number chemical potential, p, is re¬ 
lated to the baryon chemical potential by ^ = ^_b/ 3. The 
quark flelds in color, flavor and Dirac spaces are denoted 
by qia and qia = qlal^ ■ Tf are Gell-Mann matrices act¬ 
ing in flavor space. Charge conjugated quark flelds are 


denoted by = Ccf" and = q^C, where C = 17 ^ 7 ° 
is the Dirac charge conjugation matrix. The indices a, /3 
and 7 represent colors (r = 1, g = 2 and b = 3), while i, 
j and k represent flavors {u = 1, d = 2 and s = 3). Gs 
and Gd are dimensionful coupling constants that must 
be determined by experiments. 

Typically, three-flavor N.1L models use a ’t Hooft de¬ 
terminant interaction that induces a Ua( 1 ) symmetry 
breaking in the pseudoscalar isoscalar meson sector which 
can be adjusted such that the rj-rj' mass difference is de¬ 
scribed. This realization of the Uyi(l) breaking leads to 
the important consequence that the quark condensates of 
different flavor sectors get coupled. The dynamically gen¬ 
erated strange quark mass contains a contribution from 
the chiral condensates of the light flavors. There is, how¬ 
ever, another possible realization of the Ua( 1 ) symme¬ 
try breaking that does not arise on the mean field level, 
but only for the mesonic fluctuations in the pseudoscalar 
isoscalar channel. This is due to the coupling to the non- 
perturbative gluon sector via the the triangle anomaly, 
see e.g. UlliKai. This realization of the 77 - 77 ' mass 
difference gives no contribution to the quark thermody¬ 
namics at the mean held level, which we will follow in 
this paper. Up to now it is not known, which of the 
two Uu(l) breaking mechanisms is the dominant one in 
nature. In the present exploratory study of the mean 
field thermodynamics of three-flavor quark matter, we 
will take the point of view that the ’t Hooft term might 
be subdominant and can be disregarded. One possible 
way to disentangle both mechanisms is due to their dif¬ 
ferent response to chiral symmetry restoration at finite 
temperatures and densities. While in heavy-ion collisions 
only the finite temperature aspect can be systematically 
studied , the state of matter in neutron star interiors 
may be suitable to probe the Uu(l) symmetry restora¬ 
tion and its possible implications for the quark matter 
phase diagram at high densities and low temperatures. 
A comparison of the results presented in this work with 
the alternative treatment of the phase diagram of three- 
flavor quark matter including the ’t Hooft determinant 
term, see may therefore be very instructive. 

The mean-field Lagrangian is 


G — qia \_'^^dij3aP ^^G g (^qiaqj p))3ij^Sap “t“ ] qjp 

- 2GsJ2{{miV - E 

i fc,7 

^k'y ‘^G D^^H^aP'y^ijk^qi'p'j'kdjf pf'f) '^^d^ap-y^ijk^k'y- 


We define the chiral gaps 

(t>^ = -^Gs{{qiq^)), 


( 3 ) 

( 4 ) 


( 5 ) 


and the diquark gaps 

— 2GD^qia'^'yii^aP'y^ijkqjp'j) • 


( 6 ) 
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The chiral condensates contribute to the dynamical 
mass of the quarks, the constituent quark mass matrix 
in flavor space is M = diag(m° + + ^s), 

where are the current quark masses. For finite cur¬ 
rent quark masses the U(?>)l x U(?))r symmetry of the 
Lagrangian is spontaneously broken and only approxi¬ 
mately restored at high densities. 

The diquark gaps, are antisymmetric in fiavor 

and color, e.g., the condensate corresponding to A^r is 
created by green and blue down and strange quarks. Due 
to this property, the diquark gaps can be denoted with 
the fiavor indices of the interacting quarks 


Aur — Af^s, Af/q — A^ 


As 5 — Aud- 


( 7 ) 


After reformulating the mean-field lag rangian in 8 - 
component Nambu-Gorkov spinors [42, [4^ and perform¬ 
ing the functional integrals over Grassman variables ^4| 
we obtain the thermodynamic potential 




itid 


|A, 


ds 


8 Gs 


4Gz3 




(fp 


(27r)3 2 


oTrln -S' ^(iujn^p) 


T 




( 8 ) 

Here S~^{p) is the inverse propagator of the quark fields 
at four momentum p = {iuJmP), 


S {iujn,p) = 




At 


^^7 


/C7 




( 9 ) 


and u!n = (2n -I- l)7rr are the Matsubara frequencies for 
fermions. The thermodynamic potential of ultrarelativis- 
tic electrons, 

has been added to the potential, and the vacuum contri¬ 
bution. 


Do = D( 0 , 0 ) = 


blu 


4- 


i^ls 


8Gs 


(11) 


has been subtracted in order to get zero pressure in vac¬ 
uum. Using the identity Tr(ln(T>)) = ln(det(T>)) and 
evaluating the determinant (see Appendix A), we obtain 


in color, fiavor, and Nambu-Gorkov space. This result is 
in agreement with [allil- Fina lly, the Matsubara sum 
can be evaluated on closed form [43 |. 


=Aa + 2 rin(l + e-^“A)^ (14) 

T3. X 


leading to an expression for the thermodynamic potential 
on the form 


n{T,p) = 


+ + , |A,dP + 


|A,,p + |Ad.p 


8 Gs 
(fp 

Dg Dq. 


4G 


D 




(15) 


It should be noted that (14) is an even function of Aq, 
so the signs of the quasiparticle dispersion relations are 
arbitrary. In this paper, we assume that there are no 
trapped neutrinos. This approximation is valid for quark 
matter in neutron stars, after the short period of delep- 
tonization is over. 

Equations ip| . jnj, jni, and m form a consistent 
thermodynamic model of superconducting quark matter. 
The independent variables are p and T. The gaps, (j)i, 
and Aij , are variational order parameters that should be 
determined by minimization of the grand canonical ther¬ 
modynamical potential, D. Also, quark matter should be 
locally color and electric charge neutral, so at the physical 
minima of the thermodynamic potential the correspond¬ 
ing number densities should be zero 


nq = 

an 

OPQ 

(16) 

ns = 

-^ = 0 
dps ’ 

(17) 

ns = 

dps 

(18) 


The pressure, P, is related to the thermodynamic poten¬ 
tial by P = —D at the global minima of D. The quark 
density, entropy and energy density are then obtained 
as derivatives of the thermodynamical potential with re¬ 
spect to /i, T and l/T, respectively. 


III. RESULTS 


^ a—1 ^ 

( 12 ) 

The quasiparticle dispersion relations, Xaiv), are the 
eigenvalues of the Hermitian matrix, 

_ [-707 • p-fM + p i°AkjC 
[ l°CAl^ -p + ^/^M- 

( 13 ) 


The numerical solutions to be reported in this Section 
are obtained with the followingset of model parameters, 
taken from Table 5.2 of Ref. Q for vanishing ’t Hooft 
interaction, 

= 5.5 MeU , (19) 

= 112.0 MeU , (20) 

GsA^ = 2.319 , (21) 

A^ = 602.3 MeU . (22) 
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With these parameters, the following low-energy QCD 
observables can be reproduced: TO,n- = 135 MeV, uik = 
497.7 MeV, /,r = 92.4 MeV. The value of the diquark cou¬ 
pling strength Go = is considered as a free param¬ 
eter of the model. Here we present results for rj = 0.75 
(intermediate coupling) and ry = 1.0 (strong coupling). 



FIG. 1: Gaps and dynamical quark masses as a function of yr 
at T=0 for intermediate diquark coupling, rj = 0.75. 



FIG. 2: Gaps and dynamical quark masses as a function of yr 
at T=0 for strong diquark coupling, r] — 1. 


A. Quark masses and pairing gaps at zero 
temperature 

The dynamically generated quark masses and the di¬ 
quark pairing gaps are determined selfconsistently at the 
absolute minima of the thermodynamic potential, in the 
plane of temperature and quark chemical potential. This 
is done for both the strong and the intermediate diquark 
coupling strength. In Figs. 1 and 2 we show the de¬ 
pendence of masses and gaps on the quark chemical po¬ 
tential at r = 0 for r] = 0.75 and rj = 1.0, resp. A 



FIG. 3: Chemical potentials /rg and /is at T=0 for both values 
of the diquark coupling jy = 0.75 and rj — 1. All phases 
considered in this work have zero ns color charge for /is = 0, 
hence /is is omitted in the plot. 

characteristic feature of this dynamical quark model is 
that the critical quark chemical potentials where light 
and strange quark masses jump from their constituent 
mass values down to almost their current mass values do 
not coincide. With increasing chemical potential the sys¬ 
tem undergoes a sequence of two transitions: (I) vacuum 
—> two-flavor quark matter, (2) two-flavor ^ three-flavor 
quark matter. The intermediate two-flavor quark matter 
phase occurs within an interval of chemical potentials 
typical for compact star interiors. While at intermediate 
coupling the asymmetry between of up and down quark 
chemical potentials leads to a mixed NQ-2SC phase be¬ 
low temperatures of 20-30 MeV, at strong coupling the 
pure 2SC phase extends down to T=0. Simultaneously, 
the limiting chemical potentials of the two-flavor quark 
matter region are lowered by about 40 MeV. Three-flavor 
quark matter is always in the CFL phase where all quarks 
are paired. The robustness of the 2SC condensate un¬ 
der compact star constraints, with respect to changes of 
the coupling strength, as well as to a softening of the 
momentum cutoff by a formfactor, has been recently in¬ 
vestigated within a different parametrization with 
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similar trend: for rj = 0.75 and NJL formfactor the 2SC 
condensate does not occur for moderate chemical poten¬ 
tials while for 77 = 1.0 it occurs simultaneously with chiral 
symmetry restoration. Fig. 3 shows the corresponding 
dependences of the chemical potentials conjugate to elec¬ 
tric (/rg) and color (/rs) charges. All phases considered 
in this work have zero 713 color charge for 7 x 3 = 0 . 



FIG. 4: Quark-quark quasiparticle dispersion relations. For 
77 = 0.75, T = 0, and 77 = 480 MeV (left panel) there is a 
forbidden energy band above the Fermi surface. All dispersion 
relations are gapped at this point in the 77 —T plane, see Fig. 5. 
There is no forbidden energy band for the ub — sr, db — sq, 
and ur — dg — sb quasiparticles at 77 = 1, T = 84 MeV, and 
p = 500 MeV (right panel). This point in the g — T plane 
constitute a part of the gapless CFL phase of Fig. 6 . 


B. Dispersion relations and gapless phases 

In Fig. 4 we show the quasiparticle dispersion rela¬ 
tions of different excitations at two points in the phase 
diagram: (I) the CFL phase (left panel), where there is 
a finite energy gap for all dispersion relations. (II) the 
gCFL phase (right panel), where the energy spectrum is 
shifted due to the assymetry in the chemical potentials, 
such that the CFL gap is zero and (gapless) excitations 
with zero energy are possible. In the present model, this 
phenomenon occurs only at rather high temperatures, 
where the condensates are diminished by thermal fluctu¬ 
ations. 


C. Phase diagram 

The thermodynamical state of the system is character¬ 
ized by the values of the order parameters and their de¬ 
pendence on T and fi. Here we illustrate this dependency 



FIG. 5: Phase diagram of neutral three-flavor quark matter 
for intermediate diquark coupling 77 = 0.75. First-order phase 
transition boundaries are indicated by bold solid lines, while 
thin solid lines correspond to second-order phase boundaries. 
The dashed lines indicate gapless phase boundaries. The vol¬ 
ume fraction, X 2 SC, of the 2SC component of the mixed NQ- 
2SC phase is denoted with thin dotted lines, while the con¬ 
stituent strange quark mass is denoted with bold dotted lines. 

in a phase diagram. We identify the following phases: 

1. NQ: = Aas = 0; 

2 . NQ-2SC: Aud ^ 0, A„s = Ads = 0, 0<X2SC<1; 

3. 2SC: A„d ^ 0, A,, = A^s = 0; 

4. uSC: Aud ^ 0, A„s 7 ^ 0, Ads = 0; 

5. CFL: Aud ^ 0, Ada 7 ^ 0, A„s 7 ^ 0; 

and their gapless versions. The resulting phase diagrams 
for intermediate and strong coupling are given in Figs. 5 
and 6 , resp. and constitute the main result of this work, 
which is summarized in the following statements: 

1. Gapless phases occur only at high temperatures, 
above 50 MeV (intermediate coupling) or 60 MeV 
(strong coupling). 

2. CFL phases occur only at rather high chemical po¬ 
tential, well above the chiral restoration transition, 
i.e. above 464 MeV (intermediate coupling) or 426 
MeV (strong coupling). 

3. Two-flavor quark matter for intermediate coupling 
is at low temperatures (T<20-30 MeV) in a mixed 
NQ-2SC phase, at high temperatures in the pure 
2SC phase. 

4. Two-flavor quark matter for strong coupling is in 
the 2SC phase with rather high critical tempera¬ 
tures of ~ 100 MeV. 
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FIG. 6 : Phase diagram of neutral three-flavor quark matter 
for strong diquark coupling 77 = 1. Line styles as in Fig. 3 


5. The critical endpoint of first order chiral phase 
transitions is at (T, 7 r) = (44 MeV, 347 MeV) for in¬ 
termediate coupling and at (92 MeV, 305 MeV) for 
strong coupling. 

D. Quark matter equation of state 

The various phases of quark matter presented in the 
previous section have been identified by minimizing the 
thermodynamic potential, Q, in the order parameters, 
Aij and (pi. For a homogenous system, the pressure is 
P = —flmin, see Fig. 7, where the /.t-dependence of Qmm 
is shown at T = 0 for the different competing phases. 
The lowest value of flmin corresponds to the negative 
value of the physical pressure. The intersection of two 
curves corresponds to a first order phase transition. All 
other thermodynamic quantities can be obtained from 
the thermodynamic potential by derivatives. At inter¬ 
mediate coupling, we have a first order transition from 
the NQ-2SC phase to the CFL phase, whereas at strong 
coupling the first order transition is from the 2SC phase 
to the CFL phase, with a lower critical energy density. 
In Fig. 8 the equation of state for cold three-flavor quark 
matter is given on a form suitable for the investigation 
of hydrodynamic stability of gravitating compact object, 
so-called quark stars. This is the topic for the next Sub¬ 
section. 


E. Quark star configurations 

The properties of spherically symmetric, static con¬ 
figurations of dense matter can be calculated with the 
well-known Tolman-Oppenheimer-Volkoff equations for 



300 350 400 450 500 550 

p [MeV] 


FIG. 7: Minima of the thermodynamical potential for neutral 
three-flavor quark matter at T=0 as a function of the quark 
chemical potential. Note that at a given coupling 77 the state 
with the lowest flmin is attained and the physical pressure is 
P — —O ■ 

1 — ‘it'miri" 



FIG. 8 : Equation of state for three-flavor quark matter at 
T=0 with first order phase transitions. For intermediate di¬ 
quark coupling {rj = 0.75): from the mixed NQ-2SC phase to 
the CFL phase, for strong diquark coupling {tj = 1): from the 
2SC phase to the CFL phase. 

hydrostatic equilibrium of self- gravitating matter, see 
also [ 4 ^ . 

dP[r) [e^r) + P{r)][rn{r)P{r)] 
dr r[r — 2m(r)] 

Here e(r) is the energy density and P{r) the pressure at 
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distance r from the center of the star. The mass enclosed 
in a sphere with radius r is defined by 


m{r) = 47r / £{r')r'‘^dr' . 
Jo 


(24) 


These equations are solved for given central baryon 
number densities, nsir = 0 ), thereby defining a sequence 
of quark star configurations. For the generalization to 
finite temperature configurations, see Hot quark 

stars have been discussed, e.g., in In Fig. 



FIG. 9: Sequences of cold quark stars for the three-flavor 
quark matter equation of state described in the text. The ris¬ 
ing branches in the mass-central density relation (left panel) 
indicate stable compact object configurations. The mass- 
radius relations (right panel) show that the three-flavor quark 
matter described in this paper leads to very compact self¬ 
bound objects. For intermediate diquark coupling, rj — 0.75, 
stable stars consist of a mixed phase of NQ-2SC matter with a 
maximum mass of 1.21 M© (dashed line). At higher densities 
a phase transition to CFL quark matter occurs which entails 
a collapse of the star. For strong coupling, 77 = 1, the low 
density quark matter is in the 2SC phase and corresponding 
quark stars are stable up to a maximum mass of 1.326 M© 
(solid line). The phase transition to CFL quark matter entails 
an instability which at T=0 leads to a third family of stable 
stars for central densities above 9 no and a mass twin window 
of 1.301 - 1.326 Mq. 


9 we show the stable configurations of quark stars for 
the three-flavor quark matter equation of state described 
above. The obtained mass radius relations allow for very 
compact selfbound objects, with a maximum radius that 
is less than 10 km. For intermediate diquark coupling, 
rj = 0.75, stable stars consist of a NQ-2SC mixed phase 
with a maximum mass of 1.21 Mq. With incrasing den¬ 
sity, a phase transition to the CFL phase renders the se¬ 
quence unstable. For the strong diquark coupling, ?7 = 1, 
quark matter is in the 2SC phase at low densities and the 


corresponding sequence of quark stars is stable up to a 
maximum mass of 1.33 M©. The phase transition to CFL 
quark matter entails an instability that leads to a third 
family of stable stars, with masses in-between 1.30 and 
1.33 Mq. For non-accreting compact stars the baryon 



FIG. 10: Cooling an isothermal quark star configuration with 
initial mass M = 1.41 M© at temperature T = 40 MeV under 
conservation of the given baryon number N = 1.46 Nq down 
to T = 0 leads to a mass defect AM = 0.1 M© for the strong 
coupling case {r) = 1.0). Due to the twin structure at T = 0, 
two alternatives for the final state can be attained, a homo¬ 
geneous 2SC quark star or a dense 2SC-CFL quark hybrid 
star. 


number is an invariant during the cooling evolution. By 
comparing the masses of cold and hot isothermal config¬ 
urations of quark stars of equal baryon number, the max¬ 
imum mass defect (energy release due to cooling) can be 
calculated. The result for the strong diquark coupling, 
rj = 1, is shown if Fig. 10. For an initial temperature 
of 40 MeV and a given baryon number of = 1.46 Nq, 
the initial mass is M = 1.41 Mq. By cooling this object 
down to T = 0, a mass defect of AM = 0.1 Mq occurs. 
For the chosen baryon number, N = 1.46 Nq, there are 
two possible T = 0 configurations (twins). A hot star 
could thus evolve into the more compact mass-equivalent 
(twin) final state, if a fluctuation triggers the transition 
to a CFL phase in the core of the star. The structures of 
these two twin configurations are given in Fig. 11. The 
energy release of 0.1 M© is of the same order of mag¬ 
nitude as the energy release in supernova explosions and 
gamma-ray bursts. Disregarding the possible influence of 
a hadronic shell and the details regarding the heat trans¬ 
port, the cooling induced first order phase transition to 
the CFL phase could serve as a candidate process for the 
puzzling engine of these energetic phenomena 
























FIG. 11: Structure of two quark star (QS) configurations with 
M = 1.314 Mq (mass twins) for the three-fiavor quark matter 
equation of state described in the text in the case of strong 
coupling (77 = 1). The low-density twin has a radius of 9 
km and is a homogeneous 2SC quark star, the high-density 
twin is more compact with a radius of 8 km and consists of a 
CFL quark matter core with 4.65 km radius and a 2SC quark 
matter shell. 


IV. CONCLUSIONS 

We have investigated the phase diagram of three-flavor 
quark matter within an NJL model under compact star 
constraints. Local color and electric charge neutrality is 
imposed for /?—equilibrated superconducting quark mat¬ 
ter. The constituent quark masses are selfconsistently de¬ 
termined. The model refrains from adopting the’t Hooft 
determinant interaction in the mean held Lagrangian as 
a realization of the Ua( 1) symmetry breaking. Instead, 
it is assumed that the rj — r]' mass difference originates 
from an anomalous coupling of the pseudoscalar isos¬ 
inglet fluctuation to the nonperturbative gluon sector, 
which gives no contribution to the quark thermodynam¬ 
ics at the mean held level. The resulting parametrization 
of this SU/(3) NJL model results in a stronger coupling 
than NJL models with a ’t Hooft term and thus in dif¬ 
ferent phase diagrams, cf. Ref. The diquark con¬ 

densates are determined selfconsistently by minimization 
of the grand canonical thermodynamic potential. The 
various condensates are order parameters that character¬ 
ize the different phases in the plane of temperature and 
quark chemical potential. These phases are in particu¬ 
lar the NQ-2SC mixed phase, the 2SC, uSC, and CFL 
phases, as well as the corresponding gapless phases. We 
have investigated strong and intermediate diquark cou¬ 
pling strengths. It is shown that in both cases gapless 
superconducting phases does not occur at temperatures 


relevant for compact star evolution, i.e., below ~ 50 MeV. 
Three-flavor quark matter phases, e.g., the CFL phase, 
occur only at rather large chemical potential, so the ex¬ 
istence of such phases in stable compact stars is ques¬ 
tionable. The stability and stucture of isothermal quark 
star configurations are evaluated. For the strong diquark 
coupling, 2SC stars are stable up to a maximum mass of 
1.33 Mq. a second family of more compact stars (twins) 
with a CFL quark matter core and masses in-between 
1.30 and 1.33 Mq are found to be stable. For interme¬ 
diate coupling, the quark stars are composed of a mixed 
NQ-2SC phase up to a maximum mass of 1.21 Mq, where 
a phase transition to the CFL phase occurs and the con¬ 
figurations become unstable. When isothermal star con¬ 
figurations with an initial temperature of 40 MeV cools 
under conservation of baryon number, the mass defect is 
0.1 AIq for the strong diquark coupling. It is important 
to investigate the robustness of these statements, in par¬ 
ticular by inlcuding nonlocal formfactors and by going 
beyond the mean-field level by including the effects of 
a hadronic medium on the quark condensates. Finally, 
any statement concerning the occurence and stability of 
quark matter in compact stars shall include an investi¬ 
gation of the influence of a hadronic shell miziiii on 
the solutions of the equations of compact star structure. 
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APPENDIX A: DISPERSION RELATIONS 

The dispersion relations of the quasiparticles that ap¬ 
pear in the expression for the thermodynamic poten¬ 
tial ja is the eigenvalues of the Nambu-Gorkov ma¬ 
trix m- For each color and flavor combination of the 
eight component Nambu-Gorkov spinors, there is a cor¬ 
responding 8x8 entry in this matrix; for three flavors and 
three colors is a 72x72 matrix. The explicit form 
of this matrix can be represented by a table, where the 
rows and columns denote the flavor and color degrees of 
freedom 
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Qur 

Qug 

Qub 

Qdr 

Qdg 

qdb 

Qsr 

9sg 

Qsb 

Qur 

gt 
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<ll 

<ll 

1l 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Dud 

0 

0 

0 

Dus 

gt 

Hug 

0 

A 

^ug 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

-Dud 0 

0 

0 

0 

0 


0 

0 

Aub 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

~Dus 

0 

0 


0 

0 

0 

Adr 

0 

0 

0 

0 

0 

0 

~Dud 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

Adg 

0 

0 

0 

0 

Dud 

0 

0 

0 

0 

0 

0 

0 

Dds 

9db 

0 

0 

0 

0 

0 

Adb 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-Dds 0 

Qsr 

0 

0 

0 

0 

0 

0 

Asr 

0 

0 

0 

0 

~ Dus 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

Asg 

0 

0 

0 

0 

0 

0 

-Dds 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

Asb 

Dus 

0 

0 

0 

Dds 

0 

0 

0 

0 

Qur 

0 

0 

0 

0 
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0 

0 

0 

Dl 

Dur 
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0 

0 

0 

0 

0 

0 

0 

Qug 

0 

0 

0 

-Did 0 

0 

0 

0 

0 

0 

Bug 

0 

0 

0 

0 

0 

0 
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Qub 

0 

0 

0 

0 

0 

0 

-^L 

0 

0 

0 

0 

Dub 

0 

0 

0 

0 

0 

0 

Qdr 

0 

-Did 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Ddr 

0 

0 

0 

0 

0 

Qdg 

Did 

0 

0 

0 

0 

0 

0 

0 

Dl 

0 

0 

0 

0 

Bdg 

0 

0 

0 

0 

Qdb 

0 

0 

0 

0 

0 

0 

0 

-^L 0 

0 

0 

0 

0 

0 

Bdb 

0 

0 

0 

Qsr 

0 

0 

-Dl 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

D sr 

0 

0 

Qsg 

0 

0 

0 

0 

0 - 

-Dl 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Bsg 

0 

Qsb 

^us 

0 

0 

0 

Dl 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Bsb 


Each entry is a 4x4 Hermitean matrix in Dirac space. 
The diagonal submatrices are 


A, 


P + Pt.a 

0 

-M, 

0 

0 

-p + Pt,a 

0 

-M, 

-M, 

0 

-p + 

a 0 

0 

-Mi 

0 

P + Pi,a_ 


(A2) 


The eigenvalues of JAljl are the quasiparticle energies, 
Xa, that enter the thermodynamic potential 111 511 . i.e., the 
72 dispersion relations of the various quark-quark and 
antiquark-antiquark excitations. These eigenvalues can 
be calculated using a standard numerical library. How¬ 
ever, in order to reduce the computational cost, the ma¬ 
trix can be decomposed into a block-diagonal matrix by 
elementary row and column operations. 


^3,13 — 


—p — 0 Mj 0 

0 p — pj^is 0 Mj 

Mj 0 p — pj^p 0 

0 Mj 0 —p — pj^p 


(A3) 


whereas the off-diagonal blocks are given by 


D, 


0 0 0 zAij 

0 0 0 
0 i/X^p 0 0 

0 0 0 


(A4) 
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0 
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0 
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0 

0 
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0 

dL 

dL 
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dL 

0 
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0 

0 

0 

Bsb 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Hug 

0 

0 

0 

0 
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A 
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0 
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0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 
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-Dud 0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

Qsr 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-DL 
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0 

0 

0 

0 

0 

0 

nt 

Hsr 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Agr 

~ Dus 

0 

0 

0 

0 

Qub 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-DL 

Bub 

0 

0 

0 

0 

dlb 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-dL 

Bsg 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Asg 

-Dds 

Qdb 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-dL 

Bdb 


This matrix has one 24 x 24 and six 8x8 independent 
submatrices. Expressing these submatrices explicitely, 
using , the 24 x 24 matrix can be decomposed into 

two independent 12x12 submatrices by elementary row 
and column operations. Similarly, the six 8x8 matrices 


can be transformed into twelve independent 4 x 4 subma¬ 
trices. There is a two-fold degeneracy due to the Nambu- 
Gorkov basis, each matrix appears both as A 4 and 
so there are only one independent 12x12 matrix and six 
4 x 4 matrices. The 12 x 12 matrix is 


p + Pur 

0 

0 

-Mu 

0 

0 

0 

<1 

iAus 

0 

0 

0 

0 p + Pdg 

0 

0 

-Md 

0 

iAud 

0 

iAds 

0 

0 

0 

0 

0 p + Psb 

0 

0 

-Ms 

iAus 

iAds 

0 

0 

0 

0 

-Mu 

0 

0 - 

P + Pur 

0 

0 

0 

0 

0 

0 

iAud 

iAus 

0 

-Md 

0 

0 - 

P + Pdg 

0 

0 

0 

0 

'^^ud 

0 

iAds 

0 

0 


0 

0 - 

P + Psb 

0 

0 

0 

'^^us 

iAds 

0 

0 

—iAud 

iAus 

0 

0 

0 - 

P Pur 

0 

0 

Mu 

0 

0 

^Aud 

0 - 

-iAds 

0 

0 

0 

0 - 

P- Pdg 

0 

0 

Md 

0 

'^Aus — 

'iAds 

0 

0 

0 

0 

0 

0 - 

p- Psb 0 

0 

Ms 

0 

0 

0 

0 

‘^Aud 

—iAus 

Mu 

0 

0 

P f-^ur 

0 

0 

0 

0 

0 

^Aud 

0 

^Ads 

0 

Md 

0 

0 p 

— pdg 

0 

0 

0 

0 

'i'Aus 

'^Ads 

0 

0 

0 

Ms 

0 

0 p 

- Psl 


and the 4 x 4 matrices are 


M4 


P + Pi,a -Ml 0 

iAi^j -p - 0 Mj 

Mi 0 p “t" pi^ct '^Aij 

0 Mj iAij p — 


for spinor products ug — dr, ub — sr, db — sg, dr — ug, 
sr — ub, and sg — db, respectively. 

Since these matrices are Hermitean, the eigenvalues 
appear in ± pairs. Thus, in general, there are nine inde¬ 
pendent dispersion relations for quark-quark excitations 
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and nine for antiquark-antiquark excitations. The eigen¬ 
values of must be calculated numerically. The eigen¬ 
values of iFTtIi can be obtained analytically by solving for 
the roots of the quartic characteristic polynomial, 

CL^X^ CL2X‘^ uiA clq — 0, 


where 


ao = + [Mf + Mj + 

+ + MiMj + -I- 

+ MiMj + Afj - Hi^aMj - ^lj,|^Mi) , 

+ 2 (^fii^aMj — fij^fjMi p ~ ^ 

02 = ^^la + /3 - 2P" - - Mf - 2 AIj - 

O 3 — 2 {^i^a Mi,/?) ■ 

In the limit when Mi = Mj = M, which is approxi¬ 
mately valid for the ug — dr and dr — ug quasiparticles, 
the four solutions are 


A = 




A*i ,/3 


± 


g^i,a P Pj,0 


±E 


A2 


(A 9 ) 


where E = P M^. This result is in agreement 
with 0. More generally, the solutions of the quartic 
equation can be found in textbooks, see, e.g., 0 |. In 
this work the eigenvalues of the 4 x 4 matrices were calcu¬ 
lated with the exact solutions of the quartic equation and 
the eigenvalues of the 12x12 matrix were calculated with 
LAPACK. The momentum integral in m was calcu¬ 
lated with a Gaussian quadrature. The minimization of 
the thermodynamic potential was performed with conju¬ 
gate gradient methods, choosing the initial values of the 
variational parameters carefully, and then comparing the 
free energies of the various minima. The color and elec¬ 
tric charges were neutralized with a globally convergent 
Newton-Raphson method in multidimensions. 


Gapless excitations/dispersion relations are character¬ 
ized by a non-zero condensate, Aij, and a corresponding 
dispersion relation that is zero for at least one value of 
the quasiparticle momentum, i.e., the dispersion relation 
reaches the Fermi surface and there is no forbidden en¬ 
ergy band. 
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